Two-center minima in harmonic spectra from aligned polar molecules 
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We extend a model of two-center interference to include the superposition of opposite orientations 
in aligned polar molecules. We show that the position of the minimum in the harmonic spectrum 
from both aligned and oriented CO depends strongly on the relative recombination strength at 
different atoms, not just the relative phase. We reinterpret the minimum in aligned CO as an 
interference between opposite orientations, and obtain good agreement with numerical calculations. 
Inclusion of the first-order Stark effect shifts the position of the interference minimum in aligned 
CO even though aligned molecules do not posses total permanent dipoles. We explain the shift in 
terms of an extra phase that the continuum electron of oriented CO accumulates due to the Stark 
effect. 
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I. INTRODUCTION 

Two-center interference minima are an important spec- 
tral signature in high-harmonic spectra from homonu- 
clear molecules [l|, Two-center interference is often 
described within the framework of the three-step model, 
in which (i) an electron tunnels through the effective bar- 
rier formed by the molecular potential and a strong laser 
field, (ii) picks up kinetic energy in the continuum, and 
(iii) emits high-energy photons upon recombination into 
the ground state [3, 0] ■ The harmonic emission splits into 
contributions from recombination at each atomic center, 
with the possibility of destructive interference between 
different centers. The position of the minimum depends 
on the momentum of the returning continuum electron, 
and the projected internuclear distance 0, ■ 

It has recently been shown within the Lewenstein 
model (1| that two-center interference in oriented polar 
molecules is more complicated than in aligned homonu- 
clear molecules. The added complication is caused by dif- 
ferent intrinsic recombination phases at different atoms, 
which can shift the interference minimum Q. We show 
that different recombination strengths also strongly af- 
fect the minimum position, and extend the analysis to 
the case of aligned polar molecules, which have a well- 
defined alignment of the internuclear axis, but not a 
well-defined hcad-to-tail direction. Field-free molecular 
alignment is much easier to obtain experimentally than 
full orientation Q- We obtain a simple equation for the 
prediction of the minimum position by reinterpreting the 
minimum as an interference between opposite orienta- 
tions. Interference effects have also recently been consid- 
ered in homonuclear Q and in heteronuclear Q diatomic 
molecules using a different model. 

Polar molecules differ from homonuclear molecules in 



that their molecular orbitals have permanent dipole mo- 
ments that interact strongly with electric fields. This 
requires that a model of high-order harmonic generation 
(HHG) from polar molecules should include the effect of 
the laser field on the ground state, something which is 
typically ignored for atomic and homonuclear systems. 
The response of the ground state to the laser field can 
be included to lowest order by adiabatically Stark shift- 
ing the energy levels of the field-free molecular orbitals 
to the instantaneous value of the electric field 0. The 
influence of Stark shifts on strong-field ionization has re- 
cently been investigated in tunneling theory [TTl - flij ]. in 
the strong-field approximation (la ], and in calculations 
using the time-dependent Schrodinger equation in the 
single-active electron approximation 16]. 

Our present calculations show that the Stark ef- 
fect shifts the interference minimum in aligned CO by 
as much as three harmonic orders when using experi- 
mentally feasible pulse parameters even though aligned 
molecules do not posses total permanent dipoles. The 
shift is explained in terms of the phase that the contin- 
uum electron of the oriented molecule accumulates from 
the Stark effect. Our results indicate that inclusion of 
the Stark effect is important when modelling spectral fea- 
tures of HHG from polar molecules. 

The paper is organized as follows. We briefly describe 
how we calculate HHG spectra in SeaQll then present our 
model of two-center interference from aligned molecules 
in Sec. IIIII In Sec. IIVI we discuss how the Stark effect 
shifts the position of the minimum, and explain the shift 
in terms of the Stark phase accumulated by an electron 
in the continuum. We summarize our results in Sec. [V] 
Two derivations are relegated to Appendix [AJ 
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Atomic units (H = e = m c = ao = 1) are used through- 
out. 
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FIG. 1: (Color online) Parallel components of HHG spec- 
tra from CO oriented at (a) 6 = 90° and (b) 6 = 86° to 
the laser polarization axis. The smooth (red) curves are ob- 
tained by averaging the spectra over several harmonic or- 
ders. The 800 nm driving laser pulse has a peak intensity 
of 4 x 10 14 W/cm 2 . Short trajectories are selected using a 
window function, and the Stark effect is not included. 



II. CALCULATING SPECTRA 

We use CO to illustrate various aspects of HHG from 
polar molecules. For this proof-of-principle study we re- 
strict ourselves to including the highest occupied molecu- 
lar orbital (HOMO) , ignoring ionization from lower-lying 
states 
states 



J (H( 

- 20] as well as recombination into excited bound 
. The HOMO is a 3cr orbital, with vertical ion- 
ization potential I p = 14.014 eV, a permanent dipole 
fi = 1.1287 au pointing from C to O, and static polariz- 
abilities an = 3.2332 au and a± = 2.7668 au parallel and 
perpendicular to the internuclear axis (Toj . The nuclei 
are frozen to the equilibrium distance R = 1.1283 A. 
The wave function is obtained using a triple zeta va- 
lence basis set in the quantum chemistry software pack- 
age GAMESS-US 

The driving pulse is chosen with a linear polarization 
along the x axis, peak intensity 4 x 10 14 W/cm 2 , and a 
wavelength of 800 nm. The electric field has a cos 2 enve- 
lope with 10 optical cycles FWHM. The carrier envelope 
phase delay has no impact on our results for such long 
pulses, and is set to zero. 

The spectrum S a (w) of the harmonic component along 
a given direction n from a single molecule is given by 



S n (6;uj) = |A n (0;u;)| s 



(1) 



where A n (0;ui) is the projection onto n of the Fourier 
transform of the time-dependent expectation value of 
the dipole velocity operator (vdi p (#;t)) [23| . Due to 
the cylindrical symmetry of a orbitals, the molecular 
orientation is defined entirely by the angle 9 between 



FIG. 2: (Color online) Parallel components of HHG spectra 
from CO aligned at (a) 9 = 90° and (b) 6 — 86° to the laser 
polarization axis. All other parameters are the same as those 
in Fig. [U 



the laser polarization axis and the nuclear displacement 
R = R 2 — Ri. Reference [l(| explains how to calculate 
(vdi P (0;i)) using an extended Lewenstein model that in- 
cludes the time-dependent Stark shift felt by the molecu- 
lar orbitals. Short trajectories are selected with a window 
function, which suppresses trajectories that have excur- 
sion times longer than approximately two thirds of an 
optical cycle. 

Spectra for CO at two different orientations are plot- 
ted in Figure [TJ Figure EJa) shows a clear two-center 
interference minimum from CO oriented perpendicular 
to the laser polarization axis. However, as soon as the 
molecule is oriented slightly away from perpendicular in 
Fig. Gib), the minimum becomes less pronounced. Also, 
Fig. [Ub) contains even harmonics due to the broken in- 
version symmetry of the target molecule [l(| HH HH . 

The spectrum from an aligned molecule is obtained by 
adding the contribution from opposite orientations co- 
herently m, 



^aligned (#;w) = | A n (0; Uj) + A n [0 + TT] u)\ 



(2) 



The effect of alignment is seen in Fig. [2j The spectra of 
aligned and oriented CO are the same in the perpendic- 
ular geometry, but at 86° they differ substantially. The 
first observation is that the spectrum in Fig. [U[b) con- 
tains only odd harmonics, which is due to the inversion 
symmetry of the aligned target [l(| H3, The sec- 

ond observation is that the interference minimum is much 
clearer than for oriented CO. 
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FIG. 3: (Color online) Parallel components of the spectra as- 
sociated with the four trajectories of CO oriented perpendic- 
ular to the laser polarization axis. The molecular spectrum 
is the coherent sum of all four contributions, and is domi- 
nated by trajectories that start with ionization at the carbon 
atom. Short trajectories are selected, and the Stark effect is 
not included. 



III. TWO-CENTER MINIMUM IN ALIGNED 
MOLECULES 

We now present a simple model of the interference min- 
imum in Fig. [2] and show that it can be thought of ei- 
ther as an interference between recombination at differ- 
ent atoms, or as an interference between opposite orien- 
tations. 

It is a common assumption that two-center interfer- 
ence minima are completely determined by the recom- 
bination step into the ground state orbital. This as- 
sumption needs to be justified for polar molecules due 
to the different classes of electron trajectories in the con- 
tinuum [28r->30] . There are four different trajectories for 
CO, each one starting at either C or O, and recombining 
at either C or O. The contribution from each class of tra- 
jectories is plotted in Fig. [3] for CO oriented perpendic- 
ular to the laser polarization axis. Trajectories starting 
at the carbon atom completely dominate the spectrum 
due to a larger ionization dipole matrix element at car- 
bon. Removing the trajectories starting at oxygen has 
very little impact on the shape and position of the in- 
terference minimum, showing that the interference is not 
caused by the ionization step. For a detailed discussion 
of trajectories be ginning and ending on the two different 
atoms, see Rcfs. [28-30]. 

The trajectory contributions in Fig. [3] are obtained us- 
ing the extended stationary-phase approximation [28[ , in 
which the trajectories going from one atom to another ac- 
quire an extra component to the momentum. Using the 
conventional stationary-phase approximation, and thus 
neglecting this drift momentum, has a negligible effect 



on the shape and position of the interference minimum. 
Also, short trajectories were selected using a window 
function, ruling out interference between long and short 
trajectories. These observations show that continuum 
dynamics is not responsible for the interference minimum 
either. 

We now turn to the recombination step, having ruled 
out interferences in the continuum as well as in the ion- 
ization step. First the oriented molecular orbital \ip,8) 
is expanded on linear combinations of atomic orbitals 
(LCAO) \i/> n ,0) 

IV, 0) = ^c„h/> n ,0). (3) 

n 

The sum is over atomic centers. The recombination 
dipole velocity matrix element between 9) and a plane 
wave |k) then takes the form 

n 

x|c B (k|v dip |V n ,e)|. (4) 

The calculation is given in Appendix lAl The phase k-R„ 
is due to the atomic position R„, and <fi n (k.,6) is an in- 
trinsic recombination phase relating to the LCAO cen- 
tered on atom n. Equation Q shows that the harmonic 
emission from different atoms interfere, causing maximal 
destructive interference at a particular momentum k for 
a fixed angle 9. 

The key to understanding the interference in aligned 
polar molecules is to realize that the dipole velocity ma- 
trix elements of opposite orientations are closely related 

(k | v dip \ip,9 + 7r) = (k|v dip |^,f5)*. (5) 

The calculation is given in Appendix [X] Together, 
Eqs. @, ([l} and ([5]) show that the recombination dipole 
velocity matrix element of an aligned molecule is given 
by 

(k | v dip I V, fl) aligncd = Re «k | v dip 1 1/>, 6)) (6) 

= ^cos (k-R„ + 0„(k,6O) 

n 

x \cn{k\v dip \ip n ,9}\. (7) 

The position of the two-center minimum is determined 
by the value of k that minimizes the absolute value of 
Eq. (0. The simplest approach would be to ignore the re- 
combination strengths \c n (k|v d i p 9)\, and claim that 
harmonic emission from one atom cancels out the har- 
monic emission from another when 

|k|Rcos(0) + 2 (k,6») - 4>i(k,9) = (2m + 1)tt (8) 

for any integer m. In this approximation, the same result 
is obtained for oriented molecules using Eq. ^ , and was 
first presented in Ref. [6|. However, the absolute values of 
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FIG. 4: (Color online) (a) Norms of the dipole velocity ma- 
trix elements for CO oriented at 9 = 86° to k. The full 
(blue) curve is the matrix element with the entire HOMO 
(k| vdip|'!/'> 86°), the dashed (green) curve is the matrix el- 
ement with the LCAO on carbon (k|vdip|'i/'i, 86°), and the 
dash-dotted (magenta) the matrix element with the LCAO 
on oxygen (k| Vdi p |'!/'2, 86°). (b) Phases of the same matrix 
elements. The dashed and dash-dotted curves are given by 

k R„ +<t> n (k,6). 

Eqs. ((3]) and (|7|) do not generally have a minimum at the 
same location, nor can they be predicted entirely from 
recombination phases. The reason for this is that the re- 
combination strengths in Eqs. Q and (O generally have 
different |k| dependences at different atomic centers. The 
problem disappears in the case of homonuclear molecules, 
where the matrix norms are identical. 

The following example illustrates our point. Figure [4] 
shows the relevant recombination dipole velocity matrix 
elements for CO oriented at 9 = 86°. The full (blue) 
curves are the norm and the phase of the total molecular 
matrix element of CO. The dashed (green) curves are the 
norm and phase of the matrix element of the LCAO on 
carbon, and the dash-dotted (magenta) curves on oxy- 
gen. The full (blue) curve in Fig. HJa) has a minimum 
at |k| = 2.08 au. According to Eq. ([8]), the minimum 
should be located at the momentum |k| for which the 
phase difference between the carbon and oxygen matrix 
elements in Fig. E{b) is n. This gives a predicted mini- 
mum at |k| = 1.86 au, which is in poor agreement with 
the actual position. 

The interference minimum in aligned polar molecules 
can also be thought of as an interference between op- 
posite orientations. Equation ^ shows that the inter- 
ference minimum appears whenever the real part of the 
recombination matrix element of the total molecular or- 
bital (k|v < iip|'i/', 0) is zero 

Re((k|v dip |^,0» = O. (9) 

According to the solid (blue) curve in Fig. Hfb) , this hap- 



FIG. 5: (Color online) Position of the interference minimum 
(of the smoothed spectrum) from aligned CO with and with- 
out including the Stark effect in the Lewenstein model. Short 
trajectories have been selected. Also plotted are the predicted 
minima positions using a two-center interference (TCI) model 
based on Eq. @ 

pens at |k| = 2.06 au. Using the relation hcu = 0.5|k| 2 +/ p 
to convert into photon energy, this corresponds to a min- 
imum at uj — 46.3u;n. The actual position of the mini- 
mum is 47.8wn, as determined by the averaged spectrum 
in Fig. [21 For comparison, Eq. ([5]) yields a prediction of 
37.5wn- Equation (JSJ) is clearly more precise at predicting 
the minimum position than Eq. (JSJ). 

Now we are also able to understand why the interfer- 
ence minimum is stronger for aligned CO in Fig. HJb) 
than for oriented CO in Fig. [TJb). The strength of the 
interference minimum in oriented CO is determined by 
the minimal value of the matrix element norm. As il- 
lustrated in Fig. 21 this minimal value is not necessarily 
zero. For aligned CO only the real part has to be zero, 
which is obtained exactly, hence two-center interference 
minima are generally expected to be stronger for aligned 
than for oriented polar molecules. 

Equation ^ is more precise than Eq. (JSJ because no 
assumptions are made on the recombination strengths at 
different atomic centers. One might argue that it would 
be more natural to use Eq. ([7]) directly, and read off 
the minimum in the matrix element norm of the aligned 
molecule rather than focus on the matrix element phase 
of the oriented molecule. However, we show in Sec. IIVI 
that the influence of the Stark effect on the minimum 
position is more easily understood in terms of the recom- 
bination phase of the oriented molecule. 

IV. IMPORTANCE OF THE STARK EFFECT 

Equation ^ shows that the interference minimum 
from an aligned polar molecule is expected to depend on 



5 



the phase of the recombination matrix element of the ori- 
ented molecule. This begs the question whether phases 
from other steps of the HHG process influence the final 
position of the minimum. Figure [5] shows the position of 
the interference minimum in aligned CO calculated using 
an extended Lewenstein model [l0( with the pulse param- 
eters from Sec.Ql] Including the Stark effect to first order 
is seen to shift the interference minimum by as much as 
three harmonic orders. Including the second-order Stark 
shift has negligible influence on the present spectra. 

We propose that the shift of the interference minimum 
can be understood in terms of the first-order Stark shift 
that the polar HOMO experiences due to the driving 
pulse. An electron trajectory starting with ionization 
at time if , and ending with recombination at time t, will 
accumulate a first-order Stark phase [lfj given by 

$stark = / M ■ F(i")di" (10) 

= -/*cos(0)[A(t)-A(tO]. (11) 

A(t) is the vector potential of the driving field F(t) = 
— 9(A(t), which is linearly polarized along the x axis. 

The accumulated Stark phase in Eq. (Tl~Tj) can be es- 
timated using classical trajectories from the three-step 
model [!, H[. In short, we consider a single cycle of the 
electric field, and calculate the return time t and return 
energy £"kin of a classic electron being liberated at time 
t' with zero kinetic energy. The accumulated phase is 
entirely determined by t, t' and 8 through Eq. (TTT1) . and 
mapped onto photon energy using hw = E^ n + I p . The 
result is plotted as a smooth (red) curve in Fig. [5] for 
9 = 82°. The predicted Stark phase only depends on 9 
through an overall factor of cos(0). 

An alternative estimate based on the Lewenstein model 
is obtained by calculating the harmonic spectra from ori- 
ented CO with and without including the Stark effect. 
As neighbouring half-cycles experience opposite signs of 
the first-order Stark shift, we restrict the ionization step 
to times t' at which the x component of the electric field 
F(t') is positive. This singles out continuum trajectories 
that share the same sign of the Stark shift. Short tra- 
jectories are selected using a window function. The dif- 
ference in harmonic phase between the two calculations 
is plotted in Fig. [5] as an irregular (blue) curve. The 
agreement with Eq. (fTTj) is quite good, especially when 
keeping in mind the simplicity of the model, and the fact 
that there are no fitting parameters whatsoever. 

The amount by which the Stark effect shifts the in- 
terference minimum can be estimated using Eqs. @ 
and (jlip . The interference minimum appears at the mo- 
mentum |k| where the real part of the recombination ma- 
trix element is zero. For CO oriented at 9 = 86° in 
Fig. IDJb), this happens at |k| = 2.06 au, corresponding 
to uj = 46.3wn. The accumulated phase from Eq. (TTT1) is 
on the order of 0.057T, which is enough to shift the in- 
terference condition to |k| = 2.08 au, corresponding to 
u = 47.0^0- When 9 is decreased, the phase shift in- 
creases due to the cos(f5) scaling in Eq. (TTTj) . 
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FIG. 6: (Color online) Estimates of the Stark phase accumu- 
lated in the continuum for CO oriented at 8 = 82° , and using 
the pulse parameters from Sec. Mil The smooth (red) curve 
is calculated from Eq. based on classical trajectories to 
map ionization and recombination times to harmonic order. 
The irregular (blue) curve is the phase difference between har- 
monics calculated with and without including the Stark effect 
in the Lewenstein model. To avoid cancellation from neigh- 
bouring half-cycles, ionization is artificially restricted to every 
other half-cycle (see text). 

Combining Eqs. © and (ITTl) gives the two-center in- 
terference predictions plotted in Fig. [5] At perpendicular 
alignment the first-order Stark shift is exactly zero, and 
the minima coincide. The two-center model is seen to be 
off by two harmonic orders as compared with the Lewen- 
stein calculations. As the alignment angle is decreased, 
the two-center interference model adequately describes 
the shift of the minimum, down to about 8 = 80°, after 
which the Lewenstein minimum washes out and finally 
vanishes. This weakening of the interference minimum 
at smaller alignment angles is related to a slower phase 
variation of the total recombination dipole velocity ma- 
trix element, illustrated in Fig.JTJb) as a full (blue) curve 
for6» = 76°. 



V. CONCLUSIONS 

We extend a model of two-center interference to 
include the superposition of opposite orientations in 
aligned polar molecules. Our analysis shows that two- 
center interference minima can be thought of either as an 
interference between recombination at different atoms or 
as an interference between opposite orientations. Using 
the former interpretation, we get the same estimate for 
the position of the minimum as for oriented molecules [f| . 
In this approximation the minimum only depends on 
the phase difference between recombination at different 
atoms. However, calculations on CO show that the dif- 
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In the last step the intrinsic recombination phase <j> n (k, 9) 
of the matrix element c n (k | v d i p \ip n ,0) is written out 
explicitly. 

Equation ([5]) follows from the observation that the re- 
combination dipole velocity matrix element is given by 
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FIG. 7: (Color online) Same as Fig.g] but with 9 = 76°. 

ferent recombination strengths are also important in de- 
termining the position of the minimum. 

Interpreting the minimum as an interference between 
orientations, we predict the minimum position in aligned 
molecules in terms of the phase of the recombination 
dipole velocity matrix element of the oriented molecule. 
This model successfully predicts the minimum position in 
spectra calculated for aligned CO to within two harmonic 
orders for alignment angles ranging between 80 and 90 
degrees. At smaller alignment angles the observed mini- 
mum vanishes. 

Inclusion of the Stark effect shifts the observed inter- 
ference minimum even though aligned molecules do not 
posses total permanent dipoles. We explain this by the 
fact that the interference minimum is determined by the 
recombination phase of the oriented molecule, which does 
have a permanent dipole. This recombination phase is 
modified by the first-order Stark phase that the elec- 
tron picks up between ionization and recombination. Our 
model successfully reproduces the shift of the interference 
minimum due to the Stark effect. 



Appendix A 

Here we derive Eqs. ^ and ([5]). The ground state of 
the oriented molecule is denoted 9), and the corre- 
sponding wave function in the laboratory frame ip(8;r). 
Using the multi-center expansion in Eq. Q the total re- 
combination dipole velocity matrix element is 

(k | v dip \1>,0)= J ^(r)v dip V(#; r)dr (Al) 
= / ^k(r)vdi P ^2 c nipn{0; r - R„)dr 

J n 

(A2) 

= ^2 c n e lk ' r v dip V'„(6'; r - R„)dr 

(A3) 

= ^c„c ikR " / e ikr v dip V„(0;r)dr 
n J 

(A4) 



where \l){9\]&) is the Fourier transform of the ground 
state wave function. Flipping the orientation of the 
molecule is the same as adding 7r to 9. We now fix 
the coordinate system such that the laser polarization 
is along the x axis and the molecular axis confined to 
the xy plane. All molecular orbitals can then be cho- 
sen to be real, and either even or odd in z. Orbitals 
that are odd in z do not contribute to the total HHG 
spectrum 26], which leaves orbitals for which adding 
7r to 6 is the same as replacing the wave function with 
ip(9; -x, -y, z) = ip(9; -x, -y, -z). 

The recombination dipole velocity matrix element of a 
flipped molecule is then 



(k | v dip | V, 9 + tt> = kV k ^(fl + tt; k) (A8) 

= kV k ^(#; -k) (A9) 

= kV k ^*(#;k) (A10) 

= (k|v dip |V,0>*. (All) 



In the above we use the behaviour of the Fourier trans- 
form under coordinate inversion as well as a property of 
Fourier transforms of real functions. 
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